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Abstract 

The correspondence between the integrability of classical mechanical systems and 
their quantum counterparts is not a 1-1, although some close correspondencies exist. 
If a classical mechanical system is integrable with invariants that are polynomial in 
momenta one can construct a corresponding commuting set of differential operators. 
Here we discuss some 2- or 3-dimensional purely quantum integrable systems (the 
1-dimensional counterpart is the Lame equation). That is, we have an integrable po- 
tential whose amplitude is not free but rather proportional to h 2 , and in the classical 
limit the potential vanishes. Furthermore it turns out that some of these systems 
actually have N + 1 commuting differential operators, connected by a nontrivial al- 
gebraic relation. Some of them have been discussed recently by A. P. Veselov et. al. 
from the point of view of Baker- Akheizer functions. 



1 Introduction 

In classical mechanics the most common definition of integrability is that of Liouville 
integrability: Suppose we have a system with N degrees of freedom, that is, we have N 
coordinates q^ and N conjugate momenta pj, with Poisson brackets [qi,Pj] = Sij. This 
system is said to be Liouville integrable ill], if there are N functions Fk(p,q) such that 1) 
[Fi, Fj] = 0, Vi,j, and 2) the Fj are functionally independent and 3) sufficiently regular. 
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One of the F's is the Hamiltonian that gives the dynamics. Weaker and stronger types 
of integrability have also been used (e.g. partial integrability, algebraic integrability and 
super- integrability) . 

When the idea of integrability is applied to quantum mechanics the properties of com- 
mutativity and functional independence are mutually exclusive: According to a theorem 
of von Neumann 0, if a pair of self-adjoint operators A and B commutes, then there is 
a third self-adjoint operator C such that A = f(C), B = g{C). This problem is avoided 
if we concentrate on differential operators, which anyway are the most interesting from 
practical point of view. Thus, even if A = id x and B = id y commute, the operator C 
mentioned above is not a differential operator. In this paper we consider only differential 
operators and property 3) above is then automatically satisfied, but we should relax the 
independence requirement 2): we should only ignore cases where one of the commuting 
operator can be given as a polynomial of the others. 

In this article we study the existence of operators commuting with a standard type 
two-dimensional Hamiltonian (Schrodinger operator) 

H = -\h 2 (d 2 x + d 2 y ) +V{x,y). (1) 

This may be considered as arising from the classical system H = \{p 2 x + p 2 ) + V(x,y) by 
the usual correspondence p a — > —ihd a - (The classical integrability of such systems has 
been reviewed in Q.) The transition from classical to quantum integrability seems to be 
always possible, sometimes one just has to add correction terms to the Hamiltonian and to 
the operator commuting with it, but these correction terms are always 0(h 2 ) Q. Recently 
the existence of a differential operator commuting with (|]) has been discussed at length in 
|5j (along with mult i- dimensional generalizations). In that work strong assumptions were 
made on the symmetries of the potential (invariance under the Weyl group), but it turns 
out that the purely quantum integrable (PQI) systems discussed here do not have such 
symmetry properties. 

An important difference between classical and quantum integrability lies in the question 
of independence as mentioned above. Thus it will be meaningful and interesting to have 
N + 1 commuting operators, even if they must be algebraically related (sometimes this is 
called algebraic integrability). This can be already illustrated for 1-dimensional systems. 
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They are always integrable (the Hamiltonian commutes with itself) but some systems are 
more integrable that others. For example IJ the Lame operator 



C = -\h>£ i +h*P{x), (2) 

where V is the Weierstrass elliptic function, commutes with 

I=^h*^ + %ih*{±,V(x)}, (3) 

where the bracket {., .} stands for an anticommutator. The operators L and X are alge- 
braically related, 

J2 = 8 £3 _ 1^ + 1^ 3; (4) 

(where <?, are the constants characterizing V) and this relation can be used to express the 
wave-function in terms of other elliptic functions. The potential h 2 V(x) is the only one 
that has a third order commuting operator, for a fifth order operator generalizations are 
possible ||. 

The classification of commuting ordinary differential operators was studied already in 
the 1920's and more recently this problem has been studied from the point of view of 
soliton theory. The corresponding two dimensional problem (i.e., the existence of three 
commuting differential operators) was studied in ||. In the one-dimensional case the 
algebraically integrable cases are (trivially) special cases of normal (classical) integrability, 
and in || it was conjectured that this holds also in higher dimensions, but it turns out not 
to be the case. 

In Section 2 we present results from a search for a third order differential operator, 
commuting with 7i of ([!]). We find that there are indeed some PQI systems, whose potential 
is proportional to h 2 and therefore vanishes in the classical limit. In section 3 we show, that 
the integrable systems found in Section 2 also have a fourth order commuting operator, 
and that the three operators are polynomially related. In Section 4 we consider fourth 
order differential operators commuting with the Hamiltonian (|l|) with a potential that is 
a generalization of the rational of four term Calogero-Moser potential. Again some PQI 
potentials are found. 
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2 Third order commuting operator 



We will consider only those third order differential operators Z 3 whose leading order part 
has constant coefficients. The analysis of the leading terms is the same as that of the 
corresponding classical system; this follows from the Weyl correspondence and the fact 
that the leading part of the Moyal bracket (which represents the commutator) agrees with 
the Poisson bracket [Q. Since the second order operator Tt has no first order derivatives, 
we may assume that Z 3 has no second order derivatives, i.e., 

T 3 = iinf{c^l + C2%d v + cid x d 2 y + cod 3 y ) + d u d x + d w d y + d . (5) 

From the leading terms of the commutation condition [Z 3 ,7i] = one immediately derives 
a linear PDE for the potential V 

C\V X XX + (3C - 2C 2 )V XX y + (3C 3 - 2C 1 )V X yy + C 2 Vyyy = 0, (6) 

where the subscripts stand for partial derivatives. The generic solution of this equation is 

V = gi{a x x + fay) + g 2 {ot 2 x + f3 2 y) + g^a^x + (3 3 y), (7) 

where are the three solutions of cia 3 + (3co — 2c2)a 2 (3 + (3c 3 — 2ci)af3 2 + c 2 /3 3 = 0. 

(The non-generic cases need separate study.) 

The potential (0) is form invariant under coordinate rotations, which just imply changes 
in the parameters aj, This suggests writing all formulae in a rotationally covariant form. 
Furthermore, we may assume that the parameters aii, satisfy the condition 

ai + a 2 + a 3 = 0, Pi + (3 2 + p 3 = 0. (8) 

(To obtain this it is only necessary to scale the arguments of the functions g.) The sub- 
manifold (§) is invariant under parameter rotations. Instead of x, y and the corresponding 
partial derivatives let us use the quantities 

Xi = cti x + (3i y, Xi = (3i x - on y, L { = -ih(Pi d x - a { d y ), Li = -ih(ai d x + Pi d y ), 

(9) 
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which satisfy 



[Li,Xj] = -[Li,Xj] = ihDij, [Li,Xj] = [LitXj] = -zTiAy, 



(10) 



where 



Due to (|) we have D = D 12 = D 23 = D 31 = -D 2l = -D 32 = -D 13 . On the manifold © 
we have An + A 2i + A 3i = so that we have the following relations between the diagonal 
and non-diagonal elements: 



An — — (Ajj + Afcj), Aij — \{— An — Ajj + A kk ) 
where {i, j, k} is a permutation of {1, 2, 3}. Furthermore 



U + L 2 + L 3 = 0, L{A 23 + L 2 2 A 31 + L 2 3 A 12 = h 2 D 2 (d 2 x + d, 



By a direct computation we find the following expression for Z 3 



T 3 — L\A\ 3 + -^2^31 + -^3^12 
3 ( 3 



33 ( i 



k=\ 



but there is one remaining condition that can be written as 

(A 22 -A 33 )(^ 2 A n -12^i) 
+ (A 3 3-An)(^ 2 A 22 -12^ 2 ) 
+ (A n -A 22 )(^ 2 A 33 -12^ 3 ) 

- 12 [(g^'z - g[g 2 )A 12 + (g 2 g' 3 - g' 2 g 3 )A 23 + (g 3 g[ - g 3 gx) A 31 ] 



(12) 



(13) 



(14) 



(15) 



0. 

The complete solution of ([X5|) is not known, but several solutions can be found when 
we note that the Weierstrass elliptic function V satisfies the equations 



V" = 12VV' 



1 V(u) T(u)' 
1 V(v) V(v)' 
1 V(w) V(w)' 



0, where u + v + w = 0. 



(16) 
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Thus if we assume that g{ = a(P we can identify the last three terms of ( 15 ) as the above 
determinant, if 

aia 2 A 12 = a 2 a 3 A 2 3 = a 3 a x A 3X . (17) 
The first three terms of (|T5| ) imply three further equation: 

(A 22 -A 33 )(^ 2 A 11 -a 1 ) = 0, 

(A33- A n )(^ 2 A 22 -a 2 ) = 0, (18) 
(A 11 -A 22 )(^A 3 3-a 3 ) = 0. 
From ([17|) we immediately find that 

a» = A jk X, k cyclic), (19) 
so that in any case the potential will be of the form 

V cx A 23 V(a lX + M + A 31 V(a 2 x + fay) + A l2 V(a 3 x + /3 3 y). (20) 
The solutions of (0,0) can be divided into three groups: 



Case 1: All Ajj's equal. 



From (0) we find that the off-diagonal A's and therefore czj's are all equal. The equality of 



A^'s allows the parameterization aij = 2Acos(#j), $ = 2Asm(9i), an d from (|8|) it follows 
that 9 2 — Q\ + |7r, 3 = Q\ — |7r. This is the well known solution with a nonzero classical 
limit. In addition to the overall rotation there is the freedom of the amplitudes A and a. 
If we choose Q\ = it/ 6 we get the familiar form 

V3.1 = a[V(A{V3x + y))+V{A(-V3x + y))+V{-2Ay)}. (21) 
Case 2: Two Ajj's equal. 

Let us say A n = A 22 ^ A 33 . Then from (0) X = h 2 A ll /A 23 = -2ft 2 A n /A 33 and 
a i = a-2 = ^ 2 Aii,a 3 = /i 2 A 11 (2A 11 — A 33 )/A 33 . Note the overall h 2 in the amplitudes cij. 
One possible parameterization (after fixing the rotational freedom) is 

V 3 . 2 = tf A2 2 + A f [2A 2 V(Bx + Ay) + 2A 2 V(-Bx + Ay) + (B 2 - A 2 )V(-2Ay)} . (22) 



This is the solution found by Veselov et. al. ||, 10] for the algebraic special case V(x) = x 2 . 
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Case 3: All A&'s different. 

Then we must have Oj = Ti 2 Au by (pT8|), and (0) implies An + A 2 2 + A 33 = 0. Thus some 
of the parameters must be complex. One parameterizations is given by 

V 33 = h 2 A 2 [2(-l + iV3 sm(9))V(iV3Acos(9)x + A(l + iV3sm(9))y) 

+ 2(-l - iV3sm(9))V(-iV3Acos(9)x + A(l -iV3sm(6))y) 

+ 4V(-2Ay)]. (23) 

Note that even though some cti, must be complex, it is possible to get a real potential 
V, e.g., if V itself is real. This solution seems to be new. 

Note that Case 2 intersects with Case 1 at B = \^3A (a = Ah 2 A 2 ), and with Case 
3 at B = iy/3A (9 = in (J2"3"D). At these special points the system should have more 
symmetries. 

For Cases 2 and 3 the amplitudes are proportional to h 2 so they are PQI. Although the 
amplitudes are fixed, there are two nontrivial degrees of freedom as in case 1, but they are 
now all in the direction parameters 

2.1 Connection with a three dimensional system 

A Hamiltonian of the form (|1]J^) (with the term — A 12 A 23 A 31 p 2 /.D 2 added) can also be 
written three dimensional system 

h 2 

H =~Y ( A23 ^ + A31 ^ + Al2 + 9l ^ 2 - ^ + ^ ~ qi) + ^ ~ ^ (24) 

through the linear transformation 

x = -Mi + g 2 /3 2 + qM/D, 

< y = (g lQtl + q 2 a 2 + q 3 a 3 )/D, (25) 
z = (giA 13 A 12 + g 2 A 23 Ai2 + q 3 A l3 A 23 )/D 2 , 

and correspondingly 

' d x = (d 1 (3 1 A 23 + d 2 f3 2 A 13 + d 3 (3 3 A 12 )/D, 

< d y = -(d 1 a 1 A 23 + d 2 a 2 A 13 + d 3 a 3 A 12 )/D, (26) 
d z = dx + d 2 + d 3 . 

v 
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If we denote //j = Ajk, k cyclic) and use (^) the results can combined as 
H = -f (fi.dl + fi 2 d 2 2 + fi 3 d 2 3 ) + A [fi 1 V{q 2 - g 3 ) + fi 2 V{q 3 - gi) + ^V{q t - q 2 )\ , (27) 
and the three integrable cases are characterized as follows: 

1. The classical case: /ii = /i 2 = /i 3 , A free. 

2. Elliptic PQI case: /^i = /i 2 , A = —h 2 ([ii + /i 3 )/ '/ii 

3. Hyperbolic PQI case: ^ + /x 2 + ^3 = 0, A = h 2 . 

Note that for the two PQI cases the same expression A = h 2 [l — (/i! + /i 2 + /i 3 )//ii] works 
for the overall amplitude. 

In H |10j the n-dimensional elliptic case was associated to a non-Coxeter configuration 
called A n (m), similar geometric interpretation for the hyperbolic case would be interesting. 

3 Fourth order commuting operator 

In this section we consider the operator 

4 

Z 4 = (^) 4 ]Tc^- l + ... (28) 

i=0 

and now the condition of commutativity [Z^TC] = leads to the necessary condition 

< A + 2(2c - c 2 )V xxxy + 3(c 3 - cx)V xxy y + 2(c 2 - 2c 4 )K W2/ - c 3 Vy yyy = 0. (29) 

Although this is a fourth order equation it does not allow four arbitrary directions (a^, /?,), 
because of the way the q's enter in the equation. We will return to this problem later and 
first consider potential with three terms. 

3.1 Three term potential 

By direct calculation with (|7|||) one finds that the invariant can be written as 

3 

~° 2 2 -9A 2 j + i, J+ 2 + ///■ ( A/,.; . i A/,. ; . 2 + 3Aj +lii+2 )} + C(a;, ?/). (30) 
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Here the indices for A are to be taken modulo 3. The integrability condition for C can be 
written as 



where 



and 



(A 23 Li + A 31 L 2 + A 12 L 3 ) Q = (31) 

n = K 1 (A 22 -A 33 )(g^h 2 A 11 -12g' 1 g 1 ) 

+ K 2 (A 33 -A 11 )(g^fi 2 A 22 -12g' 2 g 2 ) (32) 

+ K 3 (A U - A 22 )(g 3 h 2 A 33 - 12g' 3 g 3 ) 

- 12 [(gig' 2 - 9[92)A l2 + (g 2 g' 3 - g' 2 g 3 )A 23 + (g 3 g[ - g' 3 gi)A 31 ] , 



T r 3Aj i j +1 Aj j j +2 



AijAj + i j j + 2 + 2Aj i j + i 

(The denominator of can also be written as [A 2 3l/i+A 31 Z/2+A 1 2l/3, Xi\.) The expression 
(|32] ) is the same as (^), except for the extra coefficients Ki, which however play no role if 
(|32| ) is solved by = a{P with (|T7|J18D . The integrability condition for C may have other 
solutions, but we can say at least that all solutions obtained in Section 2 also have a fourth 
order invariant. However, since 7i 2 is also a fourth order invariant commuting with TC we 
must discuss the independence of the newly found invariants. 

3.2 Relations between the invariants 

Since the system is two dimensional there can be at most two algebraically independent 
commuting quantities. Indeed on finds relations between the three operators found above, 
for example in Case 1 we have 

T 4 = -A n D 4 n\ (33) 

so that in this case the new fourth order invariant is useless. 

For the other cases the algebraic relation is less trivial and therefore the extra invariant 
provides additional information. The computations are rather extensive and we have veri- 
fied only the relation connecting the leading terms of Xj, denoted below by Ij. The results 
are as follows: 
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For case 2 when A 22 = An we have the relation 

/ 2 6 A3 3 (16A^ + 16A^A 3 3 - 20A 2 1 A 2 3 - 12A n A| + 9A 4 3 ) 
+ 12/ 2 4 / 4 A2 3 (16At 1 + 8A? X A 33 - 16A^A 33 - 2A 11 A^ 3 + 3A 4 3 ) 
+ 32/ 3 / 3 2 A 2 3 (12A 4 1 - 26A 3 1 A 33 + lSA^A 2 , - A 4 3 ) 

+ 48/ 2 / 2 A 11 A 33 (16A 3 1 - 9A n A 2 3 + 2A 3 3 ) (34) 

+ 384/ 2 / 3 2 / 4 A 11 A 33 (4A 3 1 - 9A 2 n A 33 + 6A n A 2 3 - A 3 3 ) 

+ 256/ 4 A 11 (4A 4 1 - 13Af x A33 + lSA^A 2 , - 7A U A| 3 + A 4 3 ) 

+ 64/ 3 A 2 1 (16A 2 1 - 8A n A 33 + A 2 3 ) = 0, 

which is of degree 12 in the derivatives. If An = A 33 (intersection with Case 1) this 
expression simplifies to 9A 4 1 (AnJ| + 4/4) 3 = 0, in agreement with (0). It simplifies also 
if A 33 = 4A n but this corresponds to the trivial limit B = 0. 

In case 3 we get the following complicated identity among the leading terms: 

i%a 2 -a+ l)V(a- l) 2 
-6 %\ i\{a 2 - a + l) 3 a(a - 1) (a 6 - 3a 5 + 5a 3 - 3a + 1) 
+12 it i\ i 2 (a 2 - a + l) 3 a(a - 1) (2a 6 - 6a 5 + 3a 4 + 4a 3 + 3a 2 - 6a + 2) 
+3 it i\{a 2 - a + l) 2 (3a 12 - 18a 11 + 19a 10 + 70a 9 - 140a 8 - 58a 7 + 251a 6 

- 58a 5 - 140a 4 + 70a 3 + 19a 2 - 18a + 3) 

-2 i\ i\{a 2 - a + l) 3 a(a - 1) (8a 6 - 24a 5 + 21a 4 - 2a 3 + 21a 2 - 24a + 8) 
-Ai\ % 4 (^ - a + l) 2 (2a 12 - 12a 11 - 44a 10 + 330a 9 - 425a 8 - 412a 7 + 1124a 6 

- 412a 5 - 425a 4 + 330a 3 - 44a 2 - 12a + 2) 

+2i\il(a 2 - a + l)a(a - l)(35a 12 - 210a 11 + 283a 10 + 510a 9 - 1268a 8 - 298a 7 
+ 1931a 6 - 298a 5 - 1268a 4 + 510a 3 + 283a 2 - 210a + 35) 

-18 ^^(a 2 - a + l) 2 a 2 (a - l) 2 (a 2 - a - 5)(5a 2 - 11a + 5)(5a 2 + a - 1) 

-\2i\i\i\{a 2 - a + l)a(a - l)(13a 12 - 78a 11 + 54a 10 + 445a 9 - 786a 8 - 384a 7 
+ 1485a 6 - 384a 5 - 786a 4 + 445a 3 + 54a 2 - 78a + 13) 

+3i^^(2a 18 - 18a 17 + 112a 16 - 488a 15 + 1012a 14 + 28a 13 - 3599a 12 + 4330a 11 
+ 2433a 10 - 7622a 9 + 2433a 8 + 4330a 7 - 3599a 6 + 28a 5 + 1012a 4 
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- 488a 3 + 112a 2 -18a + 2) 

+108^3 i 2 (a 2 - a + l)V(a - l) 2 (a 2 - 4a + l)(a 2 + 2a - 2)(2a 2 - 2a - 1) 

+6 u 4 it(a 2 - a + l)a(a - l)(8a 12 - 48a 11 - 120a 10 + 1040a 9 - 1323a 8 - 1476a 7 

+ 3846a 6 - 1476a 5 - 1323a 4 + 1040a 3 - 120a 2 - 48a + 8) 
-12i 4 i 2 4(a 3 - 3a + l)(a 3 - 3a 2 + l)(2a 12 - 12a 11 + 41a 10 - 95a 9 + 86a 8 + 94a 7 

- 230a 6 + 94a 5 + 86a 4 - 95a 3 + 41a 2 - 12a + 2) 

+6i 4 - 2 ) 2 (a + l) 2 (2a - l) 2 (a 2 - a + l) 2 (a - l)(a 3 - 2a 2 - a + l)(a 3 - a 2 - 2a + 1) 
+729 ifa 4 (a - l) 2 (a 2 - a + l) 2 

+4i 6 z 3 (a 2 - a + l)a(a - l)(8a 12 - 48a 11 + 195a 10 - 535a 9 + 342a 8 + 1314a 7 

- 2544a 6 + 1314a 5 + 342a 4 - 535a 3 + 195a 2 - 48a + 8) 

+2 4 z 6 (8a 18 - 72a 17 + 231a 16 - 216a 15 - 777a 14 + 3507a 13 - 6201a 12 + 3003a 11 
+ 7470a 10 - 13898a 9 + 7470a 8 + 3003a 7 - 6201a 6 + 3507a 5 - 777a 4 

- 216a 3 + 231a 2 - 72a + 8) 

-4i 3 i\{a 2 - a + l) 2 a(a - l)(a 3 - 3a + l)(a 3 - 3a 2 + 1) 
(5a 6 - 15a 5 - 3a 4 + 31a 3 - 3a 2 - 15a + 5) 
+ 4 2 (a 2 - a + l) 4 (a 3 - 2a 2 - a + l) 2 (a 3 - a 2 - 2a + l) 2 = 0, (35) 

where we have used the parameterization An = (1 — a)b, A 22 = a6, A33 = —6, and 

The conditions of Cases 2 and 3 intersect when = Ajj = —2Akk- The various 
permutations of this correspond to a = 1/2, — 1 and a = 2 and for these special values (p5|) 
becomes a square of (|3~4"D, for a = — 1 and 2 we get the relation 

i% + 12 i\ i\ - 3 %\ i\ + 36 1 2 i\ i 4 - 36 i\ + 2 i\ = (36) 

and for a = 1/2 (i.e., /xi = /x 2 = — 5/^3) 

i\ - 2\i\i\ - Yli\i\ + lUi 2 i 2 3 U - 144i 4 - 16i 3 = (37) 

(The second form is obtained from the first by 23 — * 24 — > — 2i±.) In fact one can 

show that the expression ( |35|) can be written as U 2 + (a + l) 2 (a — l/2) 2 (a — 2) 2 V = 0, 
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where U and V are some polynomials of a and the i/s. Whether fl3"o] ) factorizes also for 
some other special values is an open question. 

3.3 Four term potential 



Let us now return to equation (P9|). Since it is a fourth order linear equation is should have 
solutions of type V = J2t=i 9( a i x + P%y)- This is indeed the case, but since we really have 
only four rather than five q's at our disposal (we can use H? to eliminate c 2 , say) there 
will be a relation between a{, fy. In order to simplify subsequent computations we will use 
rotational invariance to fix ot\ = 0, and scale so that fix — 1, and otj = 1, when j ^ 1. The 
geometric relation can then be written as 

3/W4 + + (3 3 + (3 4 = 0. (38) 

In the following we will also restrict our attention only to g(x) oc 1/x 2 , i.e., 

V( x v) = — + Q2 + fl3 + fl4 (39) 
{ ' V) y 2 (x + (3 2 y) 2 + (x + (3 3 y) 2 + (x + fry)* ' ld Jj 

By direct computation we obtain the leading terms of the fourth order operator com- 
muting with the Hamiltonian P,|55|), as 



J 4 = h\d x + 4d x d 3 y p 2 p 3 p 4 - d*(fcfc + + M) 
-4h 2 [{d 2 x ,g 2 + g 3 + g 4 ] 

- {d x d y , ((/3 2 + (3 3 + p 4 ) gi + (3 2 g 2 + f3 3 g 3 + (3 A g A )} 

- {d 2 y , + fofa + M9i + Ws + ft)</2 

+ D(x,y). (40) 

This far we can proceed just with ([38|). The integrability condition for D introduces further 
conditions, and leads to the following classification: 
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Case 1 



Assume that one (3 vanishes, let us say (3± = 0. Then /3 3 = — (3 2 from (^) and we get two 
integrable potentials: First there is the well known classical one with j3 2 — ~Pz — 1 

y 2 (x + y) 2 (x — y) 2 X 2 
We note that there are two degrees of freedom, all in the amplitudes. 

Case 2 

There is another possibility with f}$ = —(3 2 arbitrary: 

^ / 1 Pl\ i2 / 1 + 1 1 \ 

= C — + ^| + h 2 ( ^ 2 + ? p- - — - — . 42 

This potential has the classical limit 1/ = Ax~ 2 + By~ 2 , with is separable and therefore has 
a quadratic second invariant (in both classical and quantum picture). The result above 
indicates that in quantum mechanics there exists a non-separable generalization for it. 
There are now two free parameters, C and (3 2) the solution presented in || corresponds to 
the choice C = Ti 2 \(m + ±) 2 , Cf3j = h 2 \(l + \) 2 

Case 3 

Next we assume that (3± ^ 0, and since due to rotational invariance we could rotate any 
of the vector to the position y, we may in fact assume that all /3's are nonzero (in other 
words that no pair of the vectors is orthogonal). From the equations it then follows that 
di = h 2 (l + Pf), which normalizes the vectors. There is one solution of this type, it is 
defined by ([38]) and 

Pi + ft + ft ~ /3 2 /3 3 - AsAt - At/?2 = 0, (43) 
which have the parameterization 

^7^t' A = 7^T TP^i' (44) 

where £ is a cubic root of unity ^ 1. This is a new result. For real a > 1 the potential 
is also real because the g% and g± terms are complex conjugates. There is only one free 
parameter, but this is probably due to the special form 1/x 2 , for example if one makes the 



same restriction in (|22|J22| ) one can scale out A. 
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4 Conclusions 



We have discussed the quantum integrability in two dimensions, and in particular the 
existence of commuting differential operators of order 3 and 4. For the third order operator 
we have considered the generic case with constant coefficients in the leading term, and found 
the classical three term Calogero-Moser system (in terms of Weierstrass functions) and two 
PQI (three term) systems, one of which was first reported in [[J. These systems also have 
a fourth order commuting operator, but only for the PQI cases is the algebraic relation 
between the three operators nontrivial. 

It is also possible to have integrable four term potentials if the commuting operator is 
of fourth order. Our analysis was here restricted to l/x 2 -type potential terms and we have 
identified two PQI examples. Whether these potentials are also algebraically integrable is 
an open question: in principle one should find for them another, still higher order invariant, 
but this seems to be a formidable task. 

With these results we have only scratched the surface of pure quantum integrability. 
What is missing in particular is a geometric characterization of the phenomena and com- 
prehensive extensions to higher orders and dimensions. The classification based on Lie 



algebra root systems |TTJ, that works so well in the classical case, does not include these 
new PQI systems. 
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